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Neutron stars in mass-transferring binaries are accreting the hydrogen 
and hehum rich matter from the surfaces of their companions. This article 
explains the physics associated with how that material eventually fuses to 
form heavier nuclei and the observations of the time dependent phenom- 
ena (such as Type I X-ray bursts) associated with the thermally unstable 
thermonuclear reactions. The majority of the article is pedagogical. Rather 
than just report from the computer output, we will (when possible) explain 
why things happen the way they do. This approach invariably comes at the 
expense of numerical accuracy but has the distinct advantage of explain- 
ing how the outcome depends on the composition of the accreting matter, 
the accretion rate and the mass, radius and thermal state of the neutron 
star. We also introduce many new analytic relations that are convenient for 
comparisons to both observations and computational results. 

After explaining nuclear burning for spherically symmetric accretion 
onto neutron stars, we discuss the possibility of asymmetric burning. In 
particular, we discuss some of the mysteries from EXOSAT observations 
of Type I X-Ray bursts and how the solution to these puzzles may lie in 
considering the lateral propagation of nuclear burning fronts around the 
star. Fully understanding this problem requires knowledge of parameters 
previously neglected such as the distribution of fresh fuel on the star, the 
magnetic field strength, and the stellar rotation. Recent iJXTi? observations 
of bursters may finally tell us some of these parameters. 



1. Introduction 

Matter accreted onto a neutron star of mass M and radius R releases 
GMrUp/ R ~ 200 MeV per nucleon, a value much larger than that re- 
leased from thermonuclear fusion (7 Mev per nucleon for fusing hydrogen 
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to helium). Our lack of precise knowledge of the global accretion rate (M) 
and M/R means that the luminosity from steady-state nuclear burning (i.e. 
at the same rate as it is accreted, L„ = EnucM, where Enuc is the nuclear 
energy released per gram) is impossible to distinguish from the accretion 
luminosity. Time dependent nuclear burning is more easily observed and is 
identified by having a time -averaged luminosity L^i. 

Hansen & Van Horn (1975) first calculated how the accreted matter 
fuses to heavier elements in steady-state, and showed that almost all mod- 
els were subject to the "thin-shell" thermal instability of Schwarzschild &; 
Harm (1965). Though they did not investigate the time dependence of the 
unstable nuclear burning, Hansen and Van Horn concluded by noting: 

"The wide range of time scales found - ranging from milliseconds to 
months, depending upon the stellar model - suggests the importance 
of a continued search for periodic or quasi-periodic phenomena over a 
similar range of time scales in the compact X-ray sources . . . ." 

Such phenomena were soon found with the discovery of recurrent, quasi- 
periodic, Type I X-ray bursts from low accretion rate (M< 10~^ Mq yr^^) 
neutron stars (Grindlay et al. 1976, Belian, Conner & Evans 1976). The 
successful association (Woosley & Taam 1976, Maraschi &: Cavaliere 1977, 
Joss 1977, Lamb & Lamb 1978) of the thermal instabilities found by Hansen 
&; Van Horn (1975) with the X-ray bursts made a nice picture of a recurrent 
cycle that consists of fuel accumulation for several hours followed by a ther- 
monuclear runaway that burns the fuel in <10 seconds. The observational 
quantity a (defined as the ratio of the time-averaged accretion luminos- 
ity to the time-averaged burst luminosity) was of the order the expected 
value (i.e. a = {GM / R) / Enuc ~ 15 — 30) for a thermonuclear origin of 
the bursts. This has proven to be one important test for attributing any 
time-dependent phenomenon to nuclear burning. As we discuss here, the 
nature of the time dependence need not always be the simple limit-cycle 
behavior of Type I X-ray bursts. 

We are focused solely on neutron stars accreting globally at rates in ex- 
cess of 10"-^'^ Mq yr~^, which is appropriate for most bright, persistent Low 
Mass X-ray binaries (in particular the "Z" and "Atoll" sources of Hasinger 
& van der Klis 1989). We begin in §2 by explaining the physics of steady- 
state nuclear burning and discussing the "thin-shell" thermal instability. 
Section 3 is a derivation of the ignition conditions and an explanation of 
the accretion rate dependent regimes of unstable burning. The energetics, 
recurrence times and durations of the bursts are discussed in §4. We do 
not review much of the pre-RXTE observations of Type I X-Ray bursts, 
as there are many reviews for that material (see Lewin, van Paradijs and 
Taam, 1995 for the most recent summary of observations of Type I X-ray 
bursts, including details of radius expansion bursts which are not discussed 
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here). However, we do summarize (in §4.2) what was learned from EXOSAT 
observations of bursters, as these pointed to the possibiUty of asymmetric 
burning. 

Section 5 opens up the discussion of asymmetric burning and presents 
what is known about the speed of propagation of burning fronts on neutron 
stars. It is often said that "X-ray Pulsars Don't Burst" and "X-ray Bursters 
don't Pulse" . Section 6 discusses this dichotomy in some detail and points 
out those accreting pulsars which might someday violate this "rule", es- 
pecially the recently discovered pulsar GRO J 1744-28 (Kouveliotou et al. 
1996, Finger et al. 1996). We close in §7 by summarizing the recent RXTE 
results on coherent oscillations during Type I X-ray bursts. 

2. Thermal Stability of Steady State Nuclear Burning 

Most neutron stars are accreting a mix of hydrogen and helium from a 
relatively unevolved companion star. The standard case we consider here 
has a hydrogen mass fraction oi X = 0.7. In this case, most of the nuclear 
energy release is from fusion of the hydrogen to helium. The very high 
temperatures reached (T > 10^ K) when this burning is unstable leads to 
production of elements near the iron group. In the steady-state burning 
case, the temperatures in the hydrogen/helium burning region are not as 
high and so the fusion to the iron group elements occurs at much larger 
depths. It is also important to keep in mind that the CNO metallicity 
(ZcNo) of the accreting material can vary by a few orders of magnitude, 
depending on what stellar population the companion is from. 

2.1. EQUATIONS AND DEFINITIONS 

Most workers use the global accretion rate, M, as the fundamental parame- 
ter for discussion of the state of nuclear burning. However, the plane-parallel 
nature of the neutron star atmosphere means that the physics of thermal 
stability and nuclear ignition actually depends on the accretion rate per 
unit area, m (Fujimoto, Hanawa Sz Miyaji 1981, hereafter FHM). This pa- 
rameter determines the local behavior on the neutron star and need not be 
the same everywhere. The local Eddington rate is 



where ctth is the Thomson scattering cross-section, rup is the proton mass 
and c is the speed of light. The ram pressure within the accreted and settling 
atmosphere is negligible, so that hydrostatic balance yields P = gy, where 
dy = —pdz is the column depth (in units of g cm~^) and g k, GM/R^ (we 
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neglect general relativistic corrections throughout this article). The local 
pressure scale height, h = P/pg, is always ^ R and so g is nearly constant 
and the atmosphere is geometrically thin. The flux is carried by radiative 
transport and, for the accretion rates of interest here, the main atmospheric 
opacity is Thomson scattering, k = Kgs = (^Th{^+X)/2mp. The flux leaving 
the plane parallel atmosphere is then 

p c d c d 

Snp dz 3k dy ' 

where a is the radiation constant. For a constant flux, this integrates to 

T' = (3) 

acg 

at pressures large relative to the photospheric value Pph ^ g/K (i.e. the 
radiative zero solution, Schwarzschild 1958). 

Steady accretion modifies the equations of particle continuity and en- 
tropy for hydrostatic settling of the atmosphere (FHM). One rewrites these 
equations in a coordinate system of fixed pressure. In this case the accreted 
matter flows through the coordinates as it is compressed by accretion of 
fresh material from above. The continuity equation for an element i (with 
number density rii) is 

^ + V.(n,.-)=Er, (4) 

where r is the sum of particle creation and destruction processes. For high 
accretion rates, there is not time for differential settling (i.e. diffusion of one 
charged species relative to another) and all elements co-move downward at 
the speed needed to satisfy mass continuity, v = m/p (Wallace, Woosley 
& Weaver 1982, Bildsten, Salpeter & Wasserman 1993). We then define a 
mass fraction Xi = pi/p = AimpUi/p where Ai is the baryon number of 
species i and expand the continuity equation to obtain 

dXi . dXi AiTUpY^r 
^ + "^^ = ^^- 

The equation for the entropy is 



r§ = -lv-F + e, (6) 

dt p 



where e is the energy release rate from nuclear burning. We write the en- 
tropy as Tds = CpT{dT/T - V^dP/P), where Vad = dluT/dhiP for an 
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adiabatic change and Cp is the specific heat at constant pressure. Then 
since the temperature can depend on both time and pressure, we find 



These are the equations that describe the hydrostatic evolution of the neu- 
tron star atmosphere and are valid while the matter accumulates on the 
neutron star. The thermonuclear instability can sometimes lead to outflow 
and hydrodynamics, for which these equations cannot be used. 

2.2. THE UPPER ATMOSPHERE BEFORE BURNING 

The neutron stars in low-mass X-ray binaries accrete via an accretion disk 
formed in the Roche lobe overflow of the stellar companion. There are 
still debates about the "final plunge" onto the neutron star surface, with 
some advocating that there is a strong enough magnetic field to control the 
final infall, while others prefer an accretion disk boundary layer. The only 
known effect on the later nuclear burning is that substantial destruction 
of all elements heavier than helium via spallation reactions can occur if 
the final plunge has an appreciable radial component (Bildsten, Salpeter &; 
Wasserman 1992). The predicted amounts of spallation are substantial and 
basically turn the accreted matter into hydrogen, helium and a mix of light 
fragments, none of them capable of sustaining a CNO cycle. 

Eventually the accreted material reaches the stellar surface, thermal- 
izes and radiates away nearly all of the ~ 200 McV nucleon"^ infall energy. 
Hence, the flux exiting the star a few scale heights beneath the thermalizing 
region is from energy released much deeper within the star, due to gravi- 
tational settling and nuclear burning. At this point, the matter is part of 
the star and is in hydrostatic balance. It is continuously compressed by the 
accretion of new material from above and eventually reaches pressures and 
temperatures adequate for thermonuclear ignition. The equations describ- 
ing the settling and compression of the the atmosphere arc heat transfer 
(equation [2]) and entropy (equation [7]). Rather than give a detailed dis- 
cussion of these "settling" solutions, let's just compare terms in the entropy 
equation in order to show what is most important. 

The compression occurs on a timcscale taccr ^ y/rh, which we compare 
to the thermal timescale set by radiative diffusion 



in order to discover whether the compression is adiabatic. Presuming that 
the atmosphere at some column depth y is carrying a flux, F, then the ratio 





tth 



AacT^ 



(8) 
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of these two timescales is 

tth CpTm 



(9) 



so that, when F S> CpTm « bkBTrn/2iimp the thermal time is much 
less than the rate of compression, where is the mean molecular weight 
of the accreting gas. This is always the case in the upper atmosphere 
(/? < 10^ g cm~^), as the energy released per proton from either steady 
burning or deeper gravitational settling is much more than fc^T in the up- 
per atmosphere. In this regime, the "trajectory" of the compressed fluid 
element is far from adiabatic as it has time to transport heat while it is 
being compressed. This simplifies the calculation, as the temperature gra- 
dient in the upper atmosphere is then set by the flux. The outer boundary 
condition plays no role.^ These settling solutions are valid until the he- 
lium burns fast enough to appreciably change either the temperature or 
the helium abundance. 

For the high accretion rates of interest here, most of the outer envelope 
is non-degenerate so P = pksT/ fiirip = gy. Wc can always correct the 
equation of state for degeneracy by reducing /x to n^ff = pkuT / Pnip < ji. 
The temperature profile is determined by the flux flowing through the at- 
mosphere, which is presumed to be a constant F = Eatrn'^n'^ where Eatm 
is the energy released per accreted gram at some large depth in the atmo- 
sphere. This parameter will be determined later and depends on whether 
the burning is in steady-state or not. Then when the flux is written as 
F = IQ'^^ Ei^m^ ergs cm~^ s~^, where Eatm = ^18 10"*^^ ergs g^^ and m = 
rhs 10^ g cm^^ s^^, the relation between the density and temperature in 
the upper atmosphere from equation (3) is 

T, = 3.32,'/»(^H!!5(1±^M)"' (10) 

where T = TglO^i^, p = p^lQ^ g cm^'^ and g = 5fi4l0^^ cm s~^.^ This gives 
a good flrst cut at how the temperature changes with depth in the accreting 
envelope above the burning region. 

As mentioned earlier, the short thermal time leads to an entropy loss 
for the compressing fluid elements. Using the above relation, wc find that 
s oc ln(T'^/^//9) oc p^/"^), and so it is clear that the entropy of the 

^The problem is much more subtle when the thermal time at some depth is comparable 
to the local compression time, y/rn. For neutron stars, this is only known to happen at 
accretion rates far in excess of the Eddington rate at depths beneath the hydrogen and 
helium burning (Bildsten & Cutler 1995, Brown & Bildsten 1997). 

use ms as the local unit for the accretion rate, which is related to the global supply 
onto the neutron star via M = ms2 X 1O"*M0 yT-^{R/10 kmf. 
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accreting material actually decreases as it is compressed to higher densities. 
This temperature profile is convectively stable and will allow us to find 
where the nuclear burning begins. The convectively stable atmosphere has 
internal buoyancy that is measured by the Brunt- Vaisala frequency (Cox 
1980) 

where h = P/pg <^ Ris the local scale height and Ti = In P/ d In p)adiab = 
5/3 is the adiabatic index. The temperature-density profile for a constant 
flux atmosphere then gives 

= ^ = ^^""P^^ (12) 
20h 20kBT ' ^ ^ 



which decreases with depth in the star since the temperature rises. The 

1/2 

frequency of local buoyant oscillations is then A'' « (50 k}l7.)gi4/Tg , com- 
parable to the inverse of the time it takes sound to travel a scale height. 



2.3. HOW HYDROGEN BURNS AT HIGH ACCRETION RATES 

For these high accretion rates, the atmospheric temperature is always in 
excess of 10'' K, so that hydrogen burns via the CNO cycle. However, at 
large densities and high temperatures, the timescale for proton captures 
becomes shorter than the subsequent f3 decay lifetimes. 

For example, consider the "'^'^N(p,7)^''^0 reaction, which is one of the 
slower steps in the cycle. The timescale for the proton capture is tcap = 
{np{av))~^ , where {av) is the thermally averaged reaction rate from Caugh- 
lan & Fowler (1988), and becomes shorter than a typical /? decay time of 
100 seconds when T > 8 x 10^ K for the relevant density of p = 10^ g cm~^. 
These temperatures are usually reached before the unstable ignition of he- 
lium takes place. The rapid proton captures makes a slightly different burn- 
ing cycle than typically encountered, which is called the hot CNO cycle 
(Hoyle & Fowler 1965) 

i2c(p, 7)^2N(p, 7)140(/3+)1^N(p, 7)15o(/3+)15n(p, a)^^C (13) 

The time to go around this catalytic hydrogen burning loop is set by the /3 
decay lifetimes of ^"^O (ti/2 = 71 s) and ^^O (ti/2 = 122 s) and is temperature 
independent and thermally stable. While this cycle is operative, all of the 
seed nuclei are locked up in ^^O and ^^O. These /3-decay limitations fix the 
hydrogen burning rate at 



eh = 5.8 X lO^^ZcNO ergs g"^ s-\ 



(14) 
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where ^CNO is the mass fraction of CNO in the accumulating matter. 

There is thus a minimum column density of fresh fuel needed on the star, 
y/j, to burn the hydrogen at the rate at which it is accreted (i.e. steady- 
state). It is found by setting ehUh = rhXEh where ~ 6.4 x 10^*^ ergs g~^ 
is the energy released from hydrogen burning to helium, giving 

Vh « 10^° g cm-^Xms ■ (15) 

V^CNO/ 

This is equivalent to the amount of matter accreted onto the star in the 
time it takes to go around the cycle enough times to consume all of the 
hydrogen, Ef^/e^ ~ (10^/Zcno)s, or about one day for typical metallicities. 
However, the matter will reach high enough temperatures within a day of 
landing on the neutron star so that the primordial helium can ignite. In this 
case helium burning occurs in a hydrogen-rich environment, which enhances 
the nuclear reaction chains and energy release (Lamb &l Lamb 1978; Taam 
& Picklum 1978, 1979; FHM; Taam 1982). 



2.4. WHEN DOES THE HELIUM START BURNING? 

Now let's ask when the helium ignites and whether or not it is thermally 
stable. The energy generation rate for helium burning to carbon is (Hansen 
& Kawaler 1994) 

esa = 5.3 X 10^1^ exp (^^) ergs g-^ s-^ (16) 

In a steady-state burning situation, the helium depletes at the depth where 
the lifetime to the nuclear reaction equals the time it takes to cross a scale 
height (Taam 1981, Fushiki & Lamb 1987) 

where E^^ = 5.84 x 10-*^^ ergs g~^ is the energy release from — C. 
Even when hydrogen is present, this condition is the appropriate one, as the 
carbon produced from helium burning increases the number of seed nuclei 
for the hot CNO cycle and allows for more rapid hydrogen consumption. 

The helium depletion condition (equation [17]) is now combined with 
the temperature profile (see equation [10]; we leave the energy released per 
accreted particle, E'lg, as a free parameter) to find the temperature at which 
the helium is burned 



ml = 68.2 



{l+XfEl^ 



exp ( . (18) 
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Solving this transcendental equation for the temperature also gives the pres- 
sure for the helium burning, and most importantly, how these quantities de- 
pend on the local accretion rate, stellar mass and radius, and abundances. In 
order to simplify the transcendental, expand the exponential about the tem- 
perature of Tg = 3.38 so that exp(-44/r8) ^ 2.22 x 10-^(r8/3.38)^^. Some 
accuracy is compromised because of this approximation, but not much. The 
resulting equation for the burning temperature is 

rj. o AO w in8T^^l/5 (-^18 + XEisf/'^^ 

Tburn = 3.43 X 10 Kmg' — -7-^^ , (19) 

= Temperature at Helium Burning Location. 

The very strong temperature sensitivity of the helium burning makes this 
temperature estimate fairly accurate. The column density at the burning 
location is 

_ 5.2 X 10'^ g cm-^ 



m 



which is a bit trickier to compare to a real calculation that finds the abun- 
dance profiles self-consistently. 

For example, a full integration from Bildsten &; Brown (1997) for a star 
with gi4 = 1.87, X = 0.7, and Y = 0.3 accreting at 7715 = 1 and burning 
in steady-state (so that Eig = 5.3) found that one-half of the accreted 
helium is burned at a location where Tg = 5.38, y = 8.3 x 10^ g cm~^ and 
p = 2x 10^ g cm~^. The matter is non-degenerate and the opacity is slightly 
less than Thomson scattering k ~ 0.4Kes ~ 0.136 cm^ g~^. Our formulae 
(equations [19] and [20]) give Tg = 5.32 and ybum = 3.3 x lO'^'g cm~^ for 
this case. As you can see the temperature is in excellent agreement. Even 
when we correct for the lower opacity, the estimated column is nearly a 
factor of 2 low relative to the numerical estimate. This is because we have 
found the place where the lifetime of an alpha particle is comparable to the 
time since arrival on the star. To completely burn the matter in steady- 
stale require, a fmo longer than this a^d .s Lghly W - lO^s/^f or 
< hour at accretion rates where the burning is stable. 

We can now check the previous statement that the helium ignites be- 
fore the hydrogen has completely burned. This is true for high accre- 
tion rates, as ybum < Vh (compare equations [15] and [20]) as long as 
7715 > 0. 016(0. Ol/ZcATo)^^^. All steady-state solutions that we find at high 
accretion rates are burning the helium in a hydrogen-rich environment. 

2.5. THERMAL STABILITY OF STEADY-STATE HELIUM BURNING 



Let's now understand whether these models are stable to a thermal pertur- 
bation. It is best to start by writing down the energy equation. Since we 
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are considering a thermal perturbation on a timescale much shorter than 
the time it takes for the matter to move a scale height, all terms oc m in 
equation (7) can be dropped, giving 

It is a good approximation to presume a constant pressure during the ther- 
mal perturbations. The simplest linear stability analysis of the steady-state 
model presumes just one zone (i.e. single pressure, Fujimoto et al. 1981) 
and so we set 

dF acT^ _ 

'd^ ~ ^ ^^^^ 

as the simple representation of the radiative cooling rate from the atmo- 
sphere. The /3-limited hydrogen burning is temperature independent, so 

that the instability must arise from helium ignition. Writing 63^ oc T'^p^, 
where v = 44/T8 — 3, a constant pressure perturbation of the equilibrium 
model (i.e. e = —€cool) is stable when 

dln/c dlnp / dlnK\ ^ 

where all the derivatives are at constant pressure. These steady-state solu- 
tions are only marginally degenerate (the density at a given pressure and 
temperature never differs from the non-degenerate guess by more than 15%) 
in the burning region and so dlnp/dlnT = —1. In addition, the opacity is 
not far from Thomson scattering, so those derivatives are zero, making the 
condition just v < A + rj for stability, clearly displaying that the thermal 
stability is a competition between nuclear heating and radiative cooling 
(which scales oc T^). We thus need T > 4.88 x 10^ K for stable helium 
burning. More sophisticated non-local linear analysis have been carried out 
by Fushiki & Lamb (1987). 

For a thin shell, the degree of degeneracy has only a small effect on the 
stability criterion and, most importantly, this instability does not require 
that the matter be degenerate. The more important condition is that it is 
thin {h <C R) so that it's temperature increases when it is heated (i.e. a 
positive specific heat, Cp > 0). It is clearly thin before burning and remains 
so even during the flash as the large gravitational well on the neutron star 
requires temperatures of order 10^^ K for h ^ R. Temperatures this high 
are never reached. This analysis provides a reasonable indicator of when we 
should expect time dependence of the burning. The strength of the "flashes" 
is sensitive to the non-linearity of the nuclear reaction rates, such as the 
weakening of the temperature dependence at higher temperatures, and the 
decreasing density as the temperature rises at fixed pressure. 
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We are lucky that it is temperature which determines stability, as this 
is the quantity that is best determined by this analytic calculation. The 
above discussion basically pointed to the condition Tju^n > 4.88 x lO^K for 
stability, which, from equation (19) turns into a condition 

.oo.inS g (Fgl4/x)V^ (t^es\'"^ .... 

m > nist = 5.83 x 10 — Ti — 77 ^ 7^77 — (24) 

cm2 s [(1 + X)Ei8f/^ \ K J ^ ^ 

for stable helium ignition. We first evaluate this for accretion of pure helium 
and then discuss in more detail the more prevalent case of accretion of 
hydrogen-rich material. In both instances, we will define an accretion rate, 
TTist, above which the helium burning is thermally stable. 

There are a few ultra-compact (Porb < 50 min) binaries (4U 1820-30, 
4U 1916-05 and 4U 1626-67, Nelson, Rappaport & Joss 1986) which are 
most likely accreting pure helium from a degenerate dwarf. The smaller 
nuclear energy released per gram for helium accretion (relative to hydrogen- 
rich accretion) means that the flux flowing through a steady-state burning 
envelope is less (Joss 1978, Joss Sz Li 1980, Bildsten 1995). Let's consider 
the steady-state helium burning case, where Y = 1, X = 0, fi = 4/3 and 
Eis = 0.58, in which case the accretion rate needed for stability is 

R g / M /lO km\ 

™,.(P>.e Helium) 1.4 xlO«j^(^) (— ) . (25) 

20 times higher than the hydrogen-rich limit we will discuss next. This is 
due to the reduced amount of energy release from helium burning and 
allows for the bright helium accreting object 4U 1820-30 (see Bildsten 
1995 for a detailed discussion) to be violently unstable at accretion rates 

2.6. STABLE HYDROGEN AND HELIUM BURNING 

Most neutron stars are accreting from relatively unevolved companions and 
so X = 0.7, Y = 0.3 and /j, = 0.6. The energy released when steadily 
burning all the hydrogen and helium to, say, carbon is Eis = 5.3. This 
gives a critical local accretion rate for stable helium burning in this mixed 
environment of 

m.,(H/He)«6.5xlO';^(— ) (— ) (^) ■ (26) 

We can also write this in terms of the global Eddington rate 



MEdd VO.3O.6 1.4M0; V 1.7 



12 



LARS BILDSTEN 



where we have inserted the typical ratio of the opacities in the burning 
region. There is no reason to believe that this calculation is accurate enough 
to say whether the critical accretion rate for burning is above or below the 
Eddington rate. What is more trustworthy are the scalings, and notice that 
equation (27) is independent of the stellar radius and weakly dependent on 
the stellar mass. Using the estimate of the critical accretion rate (cqiiation 
[27]), one finds that neutron stars brighter than Lacci- = GMMgt/R ~ 
2 X 1O^*(M/1.4M0)^/^ ergs s~^ should be stably burning the hydrogen and 
helium. Note that a more massive neutron star can have unstable nuclear 
burning at higher accretion luminosities than a less massive star. This might 
provide a possible indicator of mass for those objects where the distance is 
reasonably well known. 

Fujimoto et al. (1981) estimated Mgt (they called it Men in their Ta- 
ble 1) for 3 different neutron stars. They found Mst/{W^^MQ yr"""^) = 
15, 13, and 12 for M = 0.476, 1 and 1.41M0. Our numbers inferred from 
equation (26) are Mst/{W~^MQ yr~^) = 13, 16, and 17 for the same three 
masses and radii. These are about 20 — 30% different. The most accurate 
way to learn if the burning is thermally stable is to numerically solve the 
time dependent equations (3), (5) and (7). Ayasli & Joss (1982) and Taam, 
Woosley & Lamb (1996) performed time-dependent calculations at accre- 
tion rates comparable to the Eddington limit and found that the burning is 
still thermally unstable when M < MEdd- Ayasli & Joss (1982) also carried 
out a few super-Eddington (by a few) calculations and found thermal stabil- 
ity, in agreement with the calculations here. However, there has not been a 
thorough study (i.e. surveying the dependence of Mgt on M, R and compo- 
sition) of this problem. Indeed, little has been published on the steady-state 
burning of matter on accreting neutron stars at super-Eddington accretion 
rates. The most likely environment for this burning is on the polar caps of 
high accretion rate pulsars (see §6). 

3. Unstable Mixed Hydrogen/Helium Burning 

Since most neutron stars in LMXB's are less luminous than 10^^ ergs s^^, 
we should expect some sort of time- variable phenomena from nuclear burn- 
ing. We will find that, in one dimension, the system evolves into a limit 
cycle, accumulating fuel for some long time and then reaching a thermal 
instability which leads to observable phenomena. 

3.1. THERMONUCLEAR IGNITION CONDITIONS 

How much matter needs to be accreted before reaching an instability? In 
order to keep things simple, let's presume that the accreted matter has 
a metallicity large enough so that the flux leaving the accreting envelope 
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is just F = ChUt prior to reaching an instability, where yt is the column 
of freshly accreted matter on the star since the last thermal instability 
(Taam 1980, Wallace, Wooslcy & Weaver 1982). ^ The condition for the 
thermal instability is slightly different for these non-equilibrium solutions, 
since the entropy equation before the instability is balancing heat diffusion 
with the energy released from the hot CNO cycle and gravitational settling. 
However, neither of these energy sources can lead to a thermal instability. 
It is the helium burning which takes care of that. 

The boundary of unstable/stable helium burning in the y — T plane is 
defined by the relation 

It - 

(FHM, Hanawa & Fujimoto 1982, Fushiki & Lamb 1987), where tcooi is 
defined by equation (22). These derivatives are taken at constant pressure, 
as the instability grows faster than the pressure changes. Let's presume 
some knowledge about the answer and expand the exponential in the helium 
burning rate at Tg = 2.31 so that we get a simple relation for the ignition 
temperature on the unstable branch of the relation (equation [28]) 

1.83 X lO^K 

~ «:Vioy3/iO(^^^^y2)i/5- ^29) 

The freshly accreted matter is then linearly unstable to a thermal pertur- 
bation when T > Tign at the specified column density. This is a linear 
instability and needs no "trigger". It is just that (on the unstable branch) 
any slight upward temperature perturbation will lead to more rapid nuclear 
energy release than cooling. There is a stable branch for the ignition con- 
dition at higher temperatures (see heavy dashed line in Figure 1), where 
a positive temperature perturbation leads to more rapid cooling. The dark 
solid and dashed line in Figure 1 displays the complete ignition relation 
(equation [28]) for a neutron star of mass IAMq and i? = 10 km accreting 
matter with y = 0.3. 

Thermally unstable ignition of the accreting matter will occur if the 
settling matter reaches the unstable part of the ignition curve (dark solid) 
before burning all of the helium. Now, what is the temperature at the 
bottom of the accumulating matter? The flux in the atmosphere prior to 
the unstable helium ignition is just that generated from the stable hydrogen 
burning, F{y) = eh{yt—y), which when combined with the radiative transfer 

^This requires that thUt ^ kBT/nipm so that we can neglect the energy release from 
the compression of matter or just Zcno ^ W'^Tsrhs/ys- As we discuss in §3.3, this 
condition will be violated for interesting values of the metallicity, in which case the 
ignition conditions depend on the accretion rate, rh. 
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Figure 1. Ignition conditions on an accreting neutron star of mass M = IAMq and 
R — 10 km. The thick sohd curve is the unstable portion of hehum ignition (equation 
[29]). The atmosphere is thermally unstable to the right of this line. The heavy dashed 

line is the stable branch. The throe thin solid lines are settling solutions (from bottom 
to top) for ZcNO = 10^"^ and rhr, = 0.3, 0.75 and 3. The two thin dashed lines are for 
ZcNO = 10^'^ and rhr, = 0.15 and 0.75. The upward pointing arrow shows the direction 
in which the temperature evolves when the thermal instability occurs. 



equation gives a relation 



2ac 



(30) 



It is the intersection of this solution with the ignition condition (equation 
[29] ) that tells us the temperature 



T ~ 3 5 X 10«K^ 



and column density 



yign 



2 X 10"^ " 



g cm / 1 \ 



(31) 



(32) 



for the unstable ignition of helium in the hydrogen-rich environment. Both 
of these parameters are relatively insensitive to the value of the opacity. It 
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is important to point out that the ignition pressure, P = gyign oc g^^'^ and 
therefore the outcome of the unstable hehum flash depends on the neutron 
star mass and radius. This is in contrast to thermonuclear instabilities on 
accreting white dwarfs, where the ignition pressure is nearly independent 
of gravity, thus allowing for simple relations between the recurrence time 
and accretion rate (Shore, Livio and van den HctivcI 1994). 

The electrons are mildly degenerate at the ignition site, where the typ- 
ical density is p ~ (5 — 8) x lO^g cm~^. Putting in characteristic values 
for these quantities, namely Y = 0.3 and /x = 0.6 (we could correct for 
degeneracy to /Ug//, but the ignition column is not that sensitive to it), we 
find that the ignition conditions are 

and 

.,.«3.7X10»^( (— ) (— ) ,(34) 

where we have scaled the opacity to the values found from calculations 
which include electron conduction. This formula compares well with that 
of Fujimoto et al. (1987), who found yign oc Z'^^q. These calculations 
also tell us that the physical thickness of the accumulated layer prior to 
the instability is 10 meters, an important point for later reference. After 
discussing the different accretion rate regimes for burning, we will come 
back and use these results for predicting the nature of the resulting Type I 
X-Ray bursts. 

The three light solid lines are the envelope conditions just prior to 
unstable helium ignition for Zqno = 10^^ and local accretion rates of 
ms = 0.3, 0.75, and 3 (Bildsten & Brown 1997). These solutions should 
be viewed as upper limits to the ignition column density, as they presume 
that there is no flux at the bottom {y = yt) of the accumulating material. 
Note that, as derived above, the ignition conditions are nearly independent 
of the local accretion rate. The ignition conditions found numerically are 
yign ~ 3 X lO^g cm~2, Tign ~ 2.4 X lO^K and p 8 x 10^ gcm~^, a factor 
of 0.45 less dense than we would have guessed if the electrons were non- 
degenerate (see equations [33] and [34]). In this case, degeneracy is starting 
to play a role. The flux exiting the atmosphere is larger than just e^yign 
due to the internal gravitational energy release. 

We also show two settling solutions for a lower metallicity of Zqno = 
10~^ and 7715 = 0.15 and 0.75 (the two dashed lines). For metallicities this 
low the ignition conditions become dependent on m (see discussion in §3.3). 
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This is evident in Figure 1, where yign ~ 7 x 10^gcm~^ for rhs = 0.75 
and yign ~ 8 x lO^gcm"^ for 7715 = 0.15. Degeneracy is becoming more 
important {j-ieff ~ 0.3/u), thus changing the opacity to k ~ O.Olcm^ g~^. 
It is still the case though that our equations (34) and (33) can reasonably 
well represent what is found in the numerical simulations for the ignition. 



3.2. THE ACCRETION RATE DEPENDENCE OF UNSTABLE BURNING 

We now ask how the nature of the unstable burning depends on the local 
accretion rate. It ends up that there are two critical accretion rates that 
happen to overlap with the range of accretion rates relevant to the bright 
X-ray binaries. 

At very high rh's, the helium unstably ignites at the bottom of the 
freshly accumulated pile prior to hydrogen exhaustion, or yh > y^gn (Taam 
&; Picklum 1978, 1979, Taam 1980). This translates into a critical accretion 
rate 

[ — ^ — j iMrj • 

SO that when m > rfid the helium ignites in a hydrogen-rich environment. 
This formula agrees very well with the results of Hanawa Sz Fujimoto (1982) 
for gi4 < 2, even the scaling with the surface gravity. It also reproduces the 
critical accretion rates in Taam (1981) to within a factor of two for a range 
in metallicities and neutron star masses and radii. ^ 

For accretion rates below this value (m < rhd), the hydrogen burns 
steadily and accumulates a pure helium shell underneath it. The column 
density of hydrogen on the star is just given by equation (15), and beneath 
this a pure helium shell accumulates with time until a thermal instability 
is eventually reached. The helium ignition is always unstable. We will not 
derive the necessary conditions for the helium flash, but rather find the 
accretion rate below which the hydrogen burning shell becomes thermally 
unstable. This will happen when the temperature in the hydrogen burning 
region becomes low enough (roughly T < 8 x 10^ K, see §2.3) so that the 
proton capture reactions again determine the time to go around the CNO 
cycle loop. In this case, the temperature sensitivity of the CNO burning 
(1/ from the ^^N(p,7)^^0 reaction) leads to thermally unstable 

hydrogen burning. The temperature at the bottom of the hydrogen burning 

^Tho literature is littered with multiple designations for these critical accretion rates 
that can lead to some confusion. We are following the notation of the recent review 
article by Lewin et al. (1995). Just for the sake of clarity, note that our Md is the same 
as Msti-A.) of Fujimoto et al. (1981) and Mc2 is the same as Mst{B). For some reason 1 
and 2 are swapped in the article by Taam (1981). 
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shell (where e/^y = XEhfn) is = Snehy^ /2ac, so that T{H Shell) « 
7 X 10^ K{KX'^'m5/ZcNo)^^'^- For stable hydrogen burning, we must have 
T{H Shell) > 8 X 10^ K or m > where 

. _ 6.3 X 10^ g cm-^ s-^ / 0.04 cm^ g-^ \ /Zcivo \ .^q) 
"^""^^ X \ K ] \ 0.01 ) ' ^ ' 

For accretion rates m < mc2 the hydrogen ignition is unstable and initiates 
a flash. Our formula for mc2 reproduces the quoted numbers in Fujimoto et 
al. (1981), Taam (1981) and Wallace et al. (1982) to within a factor of two 
for a large range in mctallicities and neutron star masses and radii. This is 
reasonably good given the level of our approximations. 

The ratio of these two critical accretion rates, mci/rnc2 ~ 5 for a partic- 
ular neutron star, pointing out that the unstable ignition in a pure helium 
layer occurs only for a limited range of accretion rates. These three critical 
accretion rates that delineate at least four different states of time depen- 
dent burning are summarized in Table 1. Just for the ease of comparing 
to other results, the global rates for a 1.4M0, i? = 10 km neutron star 
accreting matter with Zcno = 10~^ and X = 0.7 are also displayed. FHM 
also outlined the three time dependent cases, calling them cases 1-3, which 
are also noted in Table 1. 



TABLE 1. Nuclear Burning Regimes at High Accretion Rates 



Range in Local Accretion Rate 

ih > ihst 

M > 2.6 X 1Q-^Mq yr-i 



ihst > rh > rhci 

2.6 X 10"® > M/(Moyr-i) > 10"'' 
(FHM Case 1) 

rhci > rh > ihc2 

10"^ > M/(M0yr-i) > 2 x 10"^° 
(FHM Case 2) 

ihc2 > rh 

2 X lO"^°M0yr-i > M 
(FHM Case 3) 



Type of Nuclear Burning 

Stable hydrogen/helium burning in a 
mixed H/Hc environment 
[Equations (24) and (26)] 

Thermally unstable helium ignition 
in a mixed H/He environment 
[Equation (35)] 

Thermally unstable pure He ignition 
after complete hydrogen burning 
[Equation (36)] 

Thermally unstable hydrogen burning 
triggers combined flash 
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3.3. HELIUM IGNITION AT HIGH ACCRETION RATES 

We initially assumed that the heat flux coming through the envelope prior 
to ignition was just that generated from the hot CNO cycle, or a value 
FcNO ^ ^hUign just prior to ignition. The typical ignition column den- 
sity found in equation (34) was then independent of the local accretion 
rate and solely a function of the metallicity and neutron star parameters. 
This approximation breaks down at high accretion rates or low metal- 
licities (possibly due to the spallation scenario of Bildsten et al. 1992) 
where other sources of energy release from the accumulating matter be- 
come important. It is easiest to see what happens by writing a flux (ac- 
tually a combination of gravitational energy release and deeper burning) 
Fdeep = Ecieepfn Siud ask what Edeep needs to be so as to compete with the 
flux in the envelope from the hot CNO cycle. This immediately tells us 
that if Edeep Z{2x IQi^ergs g-^ /m5){ZcNo/0M)^^/^^ , or only 0.2 MeV 
per accreted nucleon, then we need to reconsider the thermal state of the 
accumulating matter. Wc now go into this physics in some detail, as many 
of the interesting results from EXOSAT and RXTE are for neutron stars 
in this regime. 

Let's start by rederiving the ignition condition for F (x m. Presume 
that the flux flowing through the envelope prior to ignition is -F = Edeepin 
and leave E^eep as a free parameter. Then combining the ignition condition 
(equation [29]) with the temperature profile for an atmosphere carrying a 
constant flux F gives 

, , 1.8 X Ufg cm-2 /0.04 cnr' g-i \ ^^'^ ( 10i-\^rgs g-^ \ '^'^ 
''^-^"^^ " I ^ j l^:^ j 

(37) 

which we compare to the previous ignition column density (equation [34]) 
to find the condition on accretion rate, metallicity and E^eep that must be 
satisfied so that the ignition column depends on m. This gives 

^,.,^12lWi /aoW^y- / (^)"'- . (38) 

ms y K y V514/ V 0.01 / 

So when this condition is satisfied we should expect that the critical column 
density depends on the accretion rate in the way given in equation (37). 

Now, just how big is E^eep'^ Well, there are a few energy sources to 
consider. The simplest is the energy released as matter falls a few scale 
heights in the atmosphere (i. e. the heat lost due to the non-adiabatic grav- 
itational compression of the matter). This is roughly E^eep ~ ^ksT /ijivrip 
or 10^^ ergs g~^ at T = 3 x 10^ K. For normal metallicities, this energy 
release would only be relevant at m's so high that the burning is already 
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stable. This energy release is critical at lower metallicities and Bildsten &; 
Brown (1997) show this for a case of ZcNO = 10-^. Fi gure 1 shows a case 
for Zqno = 10~^, where this energy release is included. A much larger en- 
ergy source is the deeper burning of residual hydrogen left from a previous 
flash (AyasU & Joss 1982, Wallace et. al. 1982, Taam et al. 1993). Taam et 
al. (1993, 1996) have noted that the residual hydrogen left from a previous 
flash cannot burn until electron capture densities are reached (this is pre- 
sumably due to the slowness of proton captures on the high Z nuclei). In 
this limit a fraction Xr of the accreting hydrogen is burned at large depths, 
giving Edeep = 7.7 x 10^^ ergs g~^Xr. This energy release is so large that 
even just a few percent of unburned hydrogen can be important. 

Let's compare our simple estimates to the recent calculations of suc- 
cessive thermal instabilities by Taam et al. (1996). Their neutron star had 
M = 1.4M0, R = 9.1 km and X = 0.7. They followed the time depen- 
dent evolution for models ranging from 0.1 — 1.0 times the Eddington limit. 
None of the models went into steady-state burning, so for this case Mgt is 
clearly in excess of MecJcI- All models were time dependent and had substan- 
tial residual hydrogen, ranging at the highest m's from X,. = 0.08 — 0.21. 
Our formula gives yign ~ 10* g cm~^ over this range. Their accumulated 
columns as inferred from the time between ignitions was nearly constant at 
y !=a mtrec ~ 5 x 10^ g cm~^, and so our estimate is not so bad. 

The energy release from deeper regions leads to, on average, higher igni- 
tion temperatures (up to T = (3 — 4) x 10* K) and lower ignition pressures. 
This allows for the distinct possibility that during the accumulation prior 
to ignition, even the Hot CNO cycle is modified, as breakouts can occur at 
these temperatures via the '^^0{a, 7)"'^^Ne and ^^Ne(p, 7)^*^Na reactions (see 
Champagne & Wiescher 1992 for an overview). The lower ignition pressures 
also leads to lower peak temperatures during the bursts and typically less 
fuel burned in a particular burst (Taam et al. 1996). 

4. Time Dependent Burning and Type I X-Ray Bursts 

What really happens as time marches on and accretion continuously dumps 
matter onto the star? The simplest time-dependent behavior is a limit cycle. 
The star would accumulate fuel for a time tree = Vign/ fn until the thermal 
instability is reached, at which point the temperature rapidly rises and 
all of the fuel is burned. The first time-dependent numerical calculations 
found this type of behavior, but did not follow recurrent cycles. All that 
was found was the first burst and the ignition conditions for that burst 
were comparable to what is derived here. Before discussing in some detail 
what was learned from EXOS'jdT observations of Type I X-Ray bursts, let's 
discuss a few of the numbers that can be reliably predicted from theory. 
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4.1. BURST ENERGIES, RECURRENCE TIMES AND DURATIONS 

When in the highest accretion rate unstable regime (see Table 1), the ig- 
nition columns (see equations [34] and [37]) are in the range (5 — 20) x 
10^ g cm~^, giving recurrence times of tree ~ lO^s/ms or a few hours be- 
tween bursts for a neutron star accreting at 1O~^M0 yr~^. The burst energy 
is Ehurst ~ ^T^I^VignEnuc if all of the fuel is burned, giving 

P ^ 1 cfi'^ ( ^ \^ f Vign \ ( Enuc \ 

J^burst « X iU ergs y-:^^ j 1^108 gcm-2 ) U X 1018 ergs g-i ) ' 

(39) 

The flux from the burst is close to the Eddington limit, so that the duration 

is about the time it takes to radiate the nuclear energy at the Eddington 
limit, tijurst ~ Ejjurst/ LiEdd ^ 10 — 20 s. These three numbers, the burst 
energy, recurrence time and duration, are all in agreement with the obser- 
vations and are the "hallmarks" of our understanding that Type I X-Ray 
bursts arise from thermonuclear instabilities in the accreted matter. 

Let's briefly discuss what happens as the thermal instability develops 
into a full-scale burst. There are many papers describing the time dependent 
behavior of the burning (e.g. Taam & Picklum 1979, Ayasli &: Joss 1982, 
Wallace et al. 1982, Taam ct al. 1993, 1996) and the resulting thermonuclear 
reactions of the "rp process" which produce elements beyond iron (Wallace 
& Woosley 1981, Taam 1985, Van Wormer et al. 1994, Schatz et al. 1997). 
As emphasized by FHM, the pressure at the bottom of the accumulated pile 
is an important quantity for the evolution of the flash. This is because the 
flash occurs at fixed pressure, and eventually the matter becomes radiation 
dominated. For a typical ignition column of lO^gcm"^, the pressure is 
P = gy ^ 10^^ ergs cm~^, so ~ P gives a a maximum temperature 

Tmax ~ 1-5 X 10^ K. At this temperature, the degeneracy is lifted and the 
thermal time (from equation [8]) at the bottom of the burning material is 
about one second. In one-dimcnsional (i.e. spherically symmetric) burning, 
convection nearly always occurs before these maximum temperatures are 
reached. However, this only acts to transport the heat upwards a few scale 
heights, as the convectivc zone never reaches the photosphere (Joss 1977). 
The energy is ultimately taken out of the star via radiative diffusion. The 
physical thickness of the burning layer is if « bksT /2iJ,mpg k. 18.5 mTg, 
or about 20 meters thicker than before the thermal instability set in. 

4.2. WHAT WAS LEARNED FROM EXOSAT OBSERVATIONS 

The original information about Type I X-ray bursters came from low-Earth 
orbit satellites. These discovered the phenomena and the combination of 
theory and observation led to a firm belief that the bursts arose from ther- 
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monuclear instabilities (see Lewin et. al. 1995). However, the frequent data 
gaps and absence of long-term observations hindered studies of the accre- 
tion rate dependence of the phenomena. 

The 3.8 day orbit of EX O SAT was an excellent match for long-term 
monitoring of bursters, thus revealing the dependence of their nuclear burn- 
ing behavior on accretion rate. As noted above, in the simplest picture 
of a recurrent cycle, the time between bursts should decrease as M in- 
creases since it takes less time to accumulate the critical amount of fuel. 
Exactly the opposite behavior was observed from many low accretion rate 
(M ~ W^'^Mq yr^^) neutron stars. A particularly good example is 4U 
1705-44, where the recurrence time increased by a factor of ~ 4 when the 
accretion rate increased by a factor of « 2 (Langmeier et al. 1987, Gottwald 
et al. 1989). This is difficult to understand in terms of what we discussed 
in §3. If the star is accreting matter with Zqno = 10~^ then these accre- 
tion rates are at the boundary (M ~ Md) of unstable helium ignition in 
a hydrogen-rich environment at high M and unstable pure helium ignition 
at lower M. The expected change in burst behavior as M increases would 
then be to more energetic and more frequent bursts. This was not observed. 
The difficulties arc still present if we reduce the metallicitics via spallation. 

Other low M neutron stars showed similar behavior (see Bildsten's 
[1995] discussion of many of the "atoll" sources), each of which showed 
longer (and less periodic) recurrence times as M increased. More impor- 
tantly, the amount of fuel consumed in the bursts became much less than 
that accreted between bursts as M increased, van Paradijs, Penninx &; 
Lewin (1988) tabulated this effect for many X-ray burst sources and con- 
cluded that increasing amounts of accreting fuel are consumed in a less 
visible way than Type I X-ray bursts as M increases. The objects they 
studied have all had Eddington limited radius expansion bursts and so the 
ratio of the persistent flux to the flux during radius expansion measured the 
accretion rate in units of Eddington accretion rate. From this, one infers 
that these objects typically accrete at rates M k. (3 — 30) x 10~^° yr"-*^. 

There are six higher accretion rate "Z" sources (Sco X-1, Cyg X-2, GX 5- 
1, GX n+2, GX 340-hO, GX 349-^2, Hasingcr & van der Klis 1989), which 
are thought to be accreting near the Eddington limit of 10~^ Mq yr~^ 
(Fortner, Lamb & Miller 1989, Miller & Lamb 1992). Their nuclear burn- 
ing might be stable some of the time, as they vary in accretion rate by a 
factor of 2 — 3. None of these are known to exhibit periodic bursting be- 
havior, and only two (Cygnus X-2 and GX 17-1-2) have ever shown Type I 
X-ray bursts. The first burst found was from Cygnus X-2 (Kahn & Grind- 
lay 1984). It had a rise time < 2.56 s and a decay time « 5 — 10 s, but 
little else was certain from this observation. Kuulkers, van der Klis & van 
Paradijs (1995) found nine events from Cygnus X-2 which are reasonably 
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interpreted as Type I X-ray bursts. They all lasted about 3 seconds, had 
peak fluxes less than the persistent flux, and burst energies ~ lO^^ergs. 
However, they were so infrequent that the resulting a value was usually 
> 10^ and sometimes > 10^. In other words, these bursts were not burn- 
ing all of the accreted fuel. RXTE has seen one type I X-ray burst from 
Cygnus X-2 (Smalc 1997). Many bursts have been observed from GX 17+2 
(Tawara et al. 1984, Sztajno et al. 1986, Kuulkcrs et al. 1997) all of which 
have a > 10^ and long durations (typically 100 seconds) implying that little 
of the accreted fuel is consumed in the bursts. 

5. Beyond One-Dimensional Theories 

If Type I X-ray bursts arc the only manifestation of unstable burning, then 
frequent energetic bursts should be seen from all weakly magnetic and sub- 
Eddington neutron stars. EXOSAT did not find see this occur, nor have 
other satellites. One solution to this mystery is to assume that the accu- 
mulating matter only covers 10 % of the stellar area (say in the equatorial 
region if the matter is fed in via a thin accretion disk), thus increasing the 
local accretion rate. This might help with the Atoll sources, as the burning 
would become stable when the global rate exceeds M ~ lO~^M0yr~^. How- 
ever, it implies that the fuel stays in the equatorial region until it reaches 
ignition pressures. It is not clear that this is possible, especially since the 
accumulated matter is substantially lighter than the ashes from burning. 
This solution would not work for the Z sources, as they have shown bursts 
(though infrequent) and would be absolutely stable if the local accretion 
rate was a factor of 10 larger. One could appeal to the spherical flows and 
magnetic funneling implied by some for the QPO phenomenology (Alpar &; 
Shaham 1985, Miller & Lamb 1992) to have factors of two enhancement in 
the local accretion rate. This would allow some Z source to be stable and 
others not. 

Rather than invoking alternative accretion scenarios to explain the burst 
phenomenology, let's presume that the accumulated matter covers most of 
the star before it undergoes a thermonuclear instability and see if there is 
some physics we have previously omitted that now needs to be considered. 
The most likely candidate is the propagation of burning fronts around the 
star. 

5.1. IGNITION AND PROPAGATION OF BURNING IN BURSTS 

Once it was understood that the X-ray bursts were thermonuclear flashes, 
a few workers in the field (Joss 1978, Ruderman 1981, Shara 1982, Livio &; 
Bath 1982) pondered the possibility that the burning was not spherically 
symmetric. In this case, there would be an ignition site somewhere on the 
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star (most likely where m is highest) that would start a nuclear "fire" 
which spreads around the star by igniting the accumulated fuel. The bursts 
rise to a peak luminosity in about one second and are known (from the 
spectral evolution and energetics) to cover the whole star, thus requiring a 
lateral propagation speed of at least 10^ cm s~^. The fundamental reason 
to seriously consider this scenario is that it takes many hoTirs to accinnulate 
a thermally unstable pile of fuel, but only 10 seconds to burn it. Thus, the 
conditions must be identical to better than a part in 1000 for the local 
thermal instability to occur simultaneously over the whole star. This seems 
difficult to arrange in the accretion environment. 

One of the best quotes about what might be gained by relaxing the 
spherically symmetric assumption was from Joss (1978): 

"Our computations do not reproduce the complex and individualized 
burst structure and recurrence patterns that have been observed in some 
burst sources. It is possible that many of these complexities result from 
violations of spherical symmetry due to stellar rotation or non-spherical 
accretion. Such effects could allow a limited portion of the neutron star 
surface to undergo a flash that then ignites other regions, in a manner 
that varies from burst to burst and from one source to another. " 

Following this, a few (Fryxell & Woosley 1982, Nozakura, Ikeuchi &; Fu- 
jimoto 1984, Bildsten 1993, 1995) have estimated the lateral propagation 
speeds for burning fronts on neutron stars. The heating of the gas from 
the thermal instability occurs on a timescale of seconds, which is 10^ — 10^ 
times longer than the time it takes for sound to cross a scale height. The 
atmosphere has plenty of time to hydrostatically adjust to the changing 
temperature, thus making detonations very unlikely. So, let's begin by es- 
timating the lateral defiagration speed. 

The cold fuel ahead of the moving burning front is ignited by diffusive 
heating from the ashes. A steady burning front thus has a characteristic 
thickness, which is the distance where difi'usive heating is comparable to 
the nuclear heating I {acTj /'iK.p^ e^aY^"^ , where Tj is the undetermined 
burning temperature. The front propagates through / in the time it takes 
the nuclear burning to heat the gas to Ty, tnuc ~ CpT^ /e^a^ giving a speed 



This estimate is sensitive to the value Tj, thus requiring a more careful 
derivation than provided by (40). However, it is a good initial guess and 
gives (for P k, 10^^ erg cm~^) values of order 100 — 1000 cm s~^. Higher 
ignition pressures can increase this by a factor of ~ 10, but it seems diflacult 
(if not impossible) to reach the 10^ cm s~^ needed for bursts. 




(40) 
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Though the exact value of the deflagration speed is uncertain, there is 
an accurate lower bound that comes from considering the quenching of a 
burning front. Propagation is halted when the lateral thickness of the front, 
I, exceeds a pressure scale height, h = P/pg. In this case most of the released 
nuclear energy is radiated to vacuum rather than used to heat adjacent un- 
burned matter. Since I ~ h{tnuc/tth)^^^ , propagation requires enough fuel 
so that tnuc ~ Uh at the burning temperature, translating into a minimum 
column density (or pressure) for successful propagation (a "flammability" 
limit, which for pure helium accretion is ymin ~ 10^ g cm~^) and a mini- 
mum speed Vmin ^ h/tth ~ 100 cm s"^ (Bildsten 1995). ^ 

The propagation speed is most likely much larger when convection is 
occurring (Fryxell &; Woosley 1982). The enhanced mixing and lateral mo- 
tions will increase the propagation speeds over the laminar values. The 
convection is very efficient at transporting heat, as the time for sound to 
cross a scale height tc ~ Cg/g ~ 10~^ s is much shorter than the local 
thermal time. Steady-state convection then brings the temperature-density 
profile in the envelope very close to the adiabatic value, regardless of the 
choice of a mixing length. Crude dimensional analysis with a mixing length, 
Im gives a convective velocity 



a value comparable to what is needed for the rise time of a Type I X-ray 
burst. The fact that the turbulent scale is close to a scale height prohibits 
us from doing much more than saying that the propagation speed is com- 
parable to the convective velocity when convection occurs. 

5.2. LATERAL PROPAGATION AND OBSERVATIONS 

Now, how do we expect the propagation speeds to affect the observations? 
In the presence of convection, the high combustion speed [v ^ 10^ cm s~^) 
leads to rapid ignition (^10 s) of all connected parts of the surface which 
have enough fuel for convective combustion. This is the dominant burning 
mode at low M's, where the column densities prior to ignition are large. 
Most of the star is then convectively combustible at any one time, so that a 
recurrent thermal instability leads to quasi-periodic Type I X-ray bursting. 

However, as M increases (especially in the range above lO~^M0yr~^) 
the ignition columns are decreasing (see §3.3 and Taam et al. 1996) and 

^There is a possibility that the ignition columns become so low at high M's that the 
burning fronts are completely quenched, in which case no lateral propagation can occur. 
The resulting time dependent burning in this limit would be very difficult to observe, as 
the number of independent patches on the star would be very large. This has yet to be 
calculated for the mixed hydrogen/helium burning. 
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convection becomes less prevalent. In this regime, the slower burning speeds 
are more prevalent. This breaks the burst periodicity, reduces the amount 
of fuel available for Type I X-ray bursts, and leads to slow burning of the ac- 
cumulated matter on the neutron star. In this regime, the combustion front 
moves at a speed ~ 1000 cm s~^, so that a "ring of fire" of horizontal width 
vtth ^ 100 meters moves out from the local instability. Most importantly, 
the time for the fire to go around the star {R/v '-^ lO^s) is comparable 
to the fuel accumulation time so that only a few fires are burning at once 
(Bildsten 1995). A single fire gives a luminosity profile that roughly looks 
like 

Lfireit) ~ F{2irvt){vtth) ~ EnucVign'^'^v'^t, (42) 

where F is the local flux. When only a few fires are burning at one time, 
the nuclear burning luminosity appears as ~ hour-long flares of amplitude 
Ln ~ EnucM, or 1 — 5 % of the accretion luminosity (Bildsten 1995). 

5.3. DOES UNSTABLE NUCLEAR BURNING CAUSE THE VERY LOW 

FREQUENCY NOISE? 

Hasinger & van der Klis's (1989) analysis of ii'XOS'^T observations of the 
brightest LMXBs characterized the time dependence of the X-ray flux from 
many neutron stars. They found that excess power in the mHz-Hz range was 
often present and sometimes correlated with other source properties. This 
noise has a spectral density between flicker noise (oc 1//) and even steeper 
(oc 1//^) noise and was dubbed the "very low frequency noise" (VLFN) 
by van der Klis et al. (1987). In the time domain, it looks like ~ 1 — 5% 
luminosity variations on minute to hour long time-scales, comparable to 
what we just described might be produced at high accretion rates from the 
slow burning fires. 

The best way to decide if the VLFN has anything to do with nuclear 
burning is to see if the strength of the hour-long flares (i.e. the VLFN) 
correlates with the Type I X-ray bursts. One would not expect a correlation 
if the VLFN was from accretion fluctuations. Bildsten (1995) collated the 
EXOSAT X-ray burst data with the measured VLFN power and found a 
continuum of behavior in the brightest accreting X-ray sources. On one end 
of the continuum are the sources that always show strong VLFN and never 
show enough Type I X-ray bursts to burn all of the fuel (the Z sources), 
while on the other end are those objects which are always bursting and never 
show VLFN (the lowest accretion rate Atoll sources). The most illuminating 
objects are three (4U 1705-44, 4U 1636-53 and 4U 1735-44) which show X- 
ray bursts and VLFN simultaneously. The varying M's for these objects 
allows them to display a continuum of behavior. The VLFN is absent, or 
much reduced, when the neutron star is undergoing regular Type I X-ray 
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bursts, and increases in strength as the bursts become less fuel-consuming 
and aperiodic at higher M's. 

The correlation found by Bildsten (1995) points to the concept that 
as the accretion rate increases, there is a tradeoff of nuclear fuel between 
X-ray bursts and VLFN, favoring VLFN at high accretion rates. It also 
implies that the burning should by asymmetric on the star, allowing for 
measurements of the spin period. We discuss this in §7. 

6. The Role of a Magnetic Field 

A Type I X-ray burst has never been seen from a magnetized (B > W^^ G) 
accreting pulsar. This has led to the oft-quoted quip that "X-Ray pulsars 
don't burst" . The lack of bursts from accreting pulsars was surprising since 
they accrete at rates comparable to the X-ray bursters. Joss & Li (1980) first 
explained the lack of bursts by stabilizing the nuclear burning with a high 
local accretion rate (m > rrist)-^ The magnetic field funnels the accretion 
onto the polar cap and confines the accreted matter all the way to the 
ignition pressure. The constraint of m > rhst is easily satisfied for pulsars 
with M > W^'^Mq yr~-^, as the fractional area of the polar cap only needs 
to satisfy ^cap/4'7ri?^ < 0.01. This is well within the estimates obtained by 
either following the field lines from the magnetosphcric radius to the star 
(Lamb, Pethick, &: Pines 1973), or allowing the matter to penetrate through 
the magnetopause via a Rayleigh- Taylor instability and attach to field lines 
at smaller radii (Arons & Lea 1976; Eisner & Lamb 1977). 

However, given the range in magnetic field strengths and accretion rates, 
it seems unlikely that this could explain the absence of Type I X-ray bursts 
in all accreting pulsars. Bildsten (1995) thus suggested that, even when the 
burning is thermally unstable, a strong magnetic field might inhibit the 
rapid lateral convective motion needed for the combustion front to rapidly 
ignite the whole star (or polar cap) . The burning front would then propagate 
at the slower speed set by diffusive heat transport. In this case, the nuclear 
burning would appear as a flare with a duration of a few minutes to an hour 
(depending on the propagation speed and the lateral extent of the fuel-rich 
region). These flares would rise on a tinicscale comparable to the duration 
and not necessarily be asymmetrical in time. 

The held strength required to halt the convective (~ lO^cms"^) prop- 
agation of burning fronts is not known. Convection is potentially stabilized 
when > SnPign (Gough & Tayler 1966), which requires B>7x 10^^ G 

®Joss & Li (1980) also suggested that the reduction of the opacity in an ultra-strong 
[B > 10^^ G) field could stabilize burning for m < ihst- This is at odds with the recent 
abstract by Lamb, Miller and Taam (1996), who suggested that the reduced opacity 
would increase the critical accretion rate above which the burning becomes stable. Our 
equation (24) gave rhst oc k"^^*, and we agree with Lamb et al.'s suggestion. 
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in the helium burning regions. This is satisfied for most accreting pulsars. 
However, there are subtleties even at weaker magnetic fields, as the convec- 
tion occurs in a region with p ~ 8 x 10^ g cm~^, giving a local ram pressure 
when in steady-state convection of Pram ~ pVc ^ -P- This can most likely 
push around a magnetic field with strength B < (Svrpt;^)^/^ 5 x lO^G, but 
we do not really know how convection (and most importantly burning front 
propagation) is altered by magnetic fields in the interesting intermediate 
range pv^ < B^/Stt < P (5 x 10^ G < 5 < 7 x 10^^ G). 

We emphasize that the appearance of the instability (i.e., whether it 
looks like a Type I X-ray burst or a flare lasting a few minutes) will yield 
crucial information on the neutron star's surface magnetic field. The ac- 
creting pulsars most likely to yield time-dependent burning phenomena are 
either ones with low fields (like GRO J1744-28, Bildsten & Brown 1997) 
or very low accretion rates. Strohmayer et al. (1997) searched for Type I 
X-ray bursts from GRO J 1744-28, but could not find any. The Type I X-ray 
bursts that they did detect were from another source. Perhaps the best ob- 
jects to survey would be the faint and weakly magnetic "6-second pulsars" 
(IE 2259-^586, IE 1048.1-5937 and 4U 0142-^61, see Mereghetti & Stella 
1995). The inferred fields from their steady spin down are 

<10^2 Q and the 

global accretion rates are <1O~^^M0 yr~^. These objects could very well 
be unstably burning. The burst (or flare) recurrence time and energetics 
completely depends on whether the matter stays confined to the polar cap 
prior to ignition. Magnetohydrostatic calculations (Hameury et. al. 1983, 
Brown Sz Bildsten 1997) suggest that the matter does stay confined on the 
cap, in which case the recurrence time will be shorter than expected for 
spherical accretion and the burst energetics much smaller. Actually, the 
burst energy would just be E^urst ~ AcapEnucVign-, a value much smaller 
(potentially two to three orders of magnitude) than the typical Type I X- 
ray burst. Gotthelf & Kulkarni (1997) might have just detected such an 
event from an unknown (and very faint, Lq < 10^^ ergs s~^) object in the 
globular cluster M28. They found a very subluminous Type I X-ray burst 
with L --^ 4 x 10^^ ergs s~^ and total burst energy of 2 — 3 x 10^^ ergs, 
consistent with complete burning over a small fraction of the star. 

7. Observed Periodicities During Type I X-Ray Bursts 

There were many early (but statistically weak) indications of periodicities 
during Type I X-ray bursts that led Livio &: Bath (1982) to speculate upon 
the possibilities of rotation or surface waves as the source of periodicity. 
McDermott & Taam (1987) calculated the g-mode frequencies of the enve- 
lope during a type I X-ray burst and showed that the modes would rapidly 
change their frequency (by a factor of 2 in 10 seconds). The first convincing 
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case of a highly coherent oscihation during a Type I X-ray burst was that 
found by Schoelkopf & Kelley (1991) in Aql X-1 at 7.6 Hz. They argued 
that the periodicity was the neutron star rotation period and that the signal 
arose when the burning asymmetries were rotated in and out of view. 

Jongert Sz van der Klis (1996) searched for pulsations from 147 Type 
I X-Ray bursts detected by EXOSAT. Nothing was found, with typical 
fractional RMS amplitude limits of 2-6 % for frequencies up to (typically) 
256 Hz. The recent launch of the Rossi X-Ray Timing Explorer (RXTE) 
has dramatically changed this situation. Even though it is likely that our 
discussion of these results will be dated by the time this article appears, it 
is important to summarize what has been found so far. Table 2 contains a 
summary of what was known as of May 1997. 



TABLE 2. Rapid Periodicities During Type I X-Ray Bursts (as of May 1997) 



Object Name 


l^B (Hz) 


i^d (Hz) 


Reference 


4U 1728-34 


363 


363 


Strohmayer et. al. 1996 


KS 1731-260 


524 


260 ± 10 


Smith ot al. 1997, Wijnands & van dor Klis 


Aql X-1 


7.6,549 




Schoelkopf & Kelley 1991, Zhang et al. 1997b 


4U 1636-53 


581 


276 ± 10 


Zhang et al. 1997a, Wijnands et al. 1997 


MXB 1743-29 


589 




Strohmayer et al. 1997 



Strohmayer et al. (1996) found nearly coherent ub = 363 Hz oscillations 

during type I X-ray bursts from the atoll source 4U 1728-34. Pulsations with 
amplitudes of 2.5 — 10% were detected in six of the eight bursts analyzed at 
that time. In addition, two separate high frequency QPO's were found in the 
persistent emission. These changed with accretion rate, but maintained a 
fixed difference frequency of = 363 Hz, identical to the period seen during 
the bursts. They suggested a "beat" frequency-like model and associated 
the 363 Hz feature with the neutron star spin. 

Smith, Morgan & Bradt (1997) found ai^B = 523.92±0.05 Hz oscillation 
during the peak of a single Type I X-ray burst from KS 1731-260. This burst 
had radius expansion and the periodicity was only found at the end of the 
contraction phase. The oscillation was very coherent (Q > 900) and had a 
pulse fraction of 6%. No harmonics were detected. Wijnands & van der Klis 
(1997) found kHz QPO's from this object at a fixed difference frequency 
frf = 260 lb 10 Hz, consistent with one-half the frequency seen during the 
bursts. Another object, 4U 1636-53, has a difference frequency in the kHz 
QPO's of i^d = 276 ± 10 Hz (Wijnands et al. 1997) and a burst frequency 
of ub = 581 Hz (Zhang et al. 1997a), not quite consistent with twice the 
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difference frequency. Strohmayer et al. detected 589 Hz oscillations during 
three Type I X-ray bursts from an object (most likely MXB 1743-29) near 
the galactic center. These bursts were temporally separated by just over a 
month and the oscillations were detected in the few seconds after the burst 
peak. 

The frequency during the bursts from both MXB 1743-29 and 4U 1728- 
34 increased by a few Hz as the neutron star atmosphere cools. There is 
no torque large enough to change the neutron star spin this rapidly, so 
Strohmayer et al. (1997) have argued that the observed periodicity is the 
rotation rate of the burning shell. They noted that the slight (Ar ^ R) 
hydrostatic expansion during the burning can explain the observations if the 
shell conserves angular momentum. In that case, the change in thickness 
(estimated from §'s 3.1 and 4.1 to be Ar « 20 meters) will lead to a 
frequency change of the burning material by an amount Ai^ !=a i'(Ar/R) ?a 
2 X 10~^v, or a change of 1 Hz for a 500 Hz rotation. This is close to what is 
observed. It is presumed that the neutron star spin frequency is the higher 
value and is, of course, unchanging throughout the burst. 

However, the observed frequency evolution is always from low to high 
during the cooling phase of the burst. Why is this? It is most likely because 
the atmospheric scale height is decreasing during the cooling phase. The 
atmosphere expanded and spun-down at the onset of the thermonuclear 
instability. However, this spin-down is difficult to observe, as the radiative 
layer (see discussion in §4.1) delays the information about the burst, so that 
by the time the observer sees it, the expansion and spin-down has already 
occurred. There might be some bursts where the spin-down evolution can be 
detected. Such a detection would obviously make us much more confident 
of Strohmayer et al's (1997) interpretation. 

Their conjecture implies that the burning matter wraps around the star 
a few times during the instability, as the relative shear is a few Hz for 
many seconds. Shear layers are typically unstable to the Kelvin-Helmholtz 
instability. However, the buoyancy due to the mean molecular weight con- 
trast or even thermal (i.e. the buoyancy in equation [12]) can easily stabi- 
lize this shear (at least for a dynamical time) as the Richardson number 
Ri = N"^ /{dvrot/dz)'^ > 1/4 (Fujimoto 1988). The shear might thus persist 
for the few seconds required. Further theoretical studies need to be car- 
ried out to fully understand the repercussions of this result. Problems that 
immediately come to mind are the role of any magnetic field and more im- 
portantly, the possibility of longer timescale rotational instabilities (such as 
the baroclinic instability) . Basically, one must rule out any mechanism that 
can transport angular momentum to the shell in less than ten seconds. That 
is a very long time on the neutron star, about 10^ sound crossing times! 
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